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A general and computable criterion for fc- (in) separability in continuous multipartite quantum sys- 
tems is presented. The criterion can be experimentally implemented with a finite and comparatively 
low number of local observables. We discuss in detail how the detection quality can be optimised. 



I. INTRODUCTION 



Quantum entanglement has been studied quite intensely over the last few decades, resulting in a rather wide 
understanding of simple entangled systems (i.e. bipartite systems of low dimensions, especially two qubit systems; 
for an overview, see e.g. Ref. However, still many puzzling features and open questions are revealed in more 

general systems, such as multipartite systems or systems of high (or, in particular, infinite) dimensions. 
In multipartite entangled systems (which are of grave importance to technological applications of quantum informa- 
tion theory, such as quantum secret sharing or quantum computation Q), complications arise (among others) due 
to the multiple different forms in which a multipartite state can be entangled (see e.g. [1, In particular, while a 
bipartite state is either entangled or separable, a multipartite state can be partially entangled (as investigated e.g. 
in Refs [1,0]), as opposed to genuinely multipartite entangled (see e.g. [sl-fT^). 

In infinite-dimensional (continuous) systems, problems arise because certain notions of finite-dimensional systems are 
no longer met (see e.g. [l^, Nevertheless, many concepts have been generalised from the finite-dimensional case 
to the continuous one during the last decade, most noteworthy the PPT-criterion [l5| and the scheme of quantum 
teleportation [16j] (which has also already been experimentally verified il7i|). 

Rather seldomly are systems considered which contain both these sets of difficulties (e.g. [3 13 1 although 
multipartite continuous quantum systems do hold the possibility for certain applications, such as certain kinds of 
teleportation networks ^20,] . 

In this letter, we use a general framework (which was introduced in Refs. [Ill, [2^ for finite-dimensional systems) to 
formulate a criterion for partial separability (fc-separability) of arbitrary states of a continuous variable multipartite 
system, thus proving implicitly, that the framework also works perfectly well in infinite-dimensional systems. 
The article is organised as follows. In section |TT1 the basic definitions and terminology will be reviewed, such that 
in section IIIII we can present our criterion for continuous variable fc-separability, which is the main result of this 
letter. In section HVl a guideline to the application of the criterion will be given (such that its detection power can 
be optimally used). The criterion and its application are then demonstrated in two exemplary cases in section [V] 
Finally, in section IVTl we show how the criterion can also be implemented experimentally. 
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II. BASIC DEFINITIONS 



In order to formulate our result, we firstly need to define the concept of fc-separability for continuous variable 
systems. A general n-partite pure quantum state can be written in the form 

/oo 
d"x *(a;i, X2, • • • , a;„) \xi) «) jxa) • • • |a;„) . (I) 
-oo 

It is called fc-separable (fc < n) iff its distribution function factorises into k factors, i.e. 

*(xi,X2,--- ,xn) = *i({a;»J) • *2({a;«J) (2) 

where {xi } denotes a set of coordinates, i.e. corresponds to one or several particles. That is, a pure state is called 
fc-separable iff there is a fc-partition ({a^iilKa^isll ' ' ' K^i*;}) with respect to which it is separable. If fc = n then the 
state is called fully separable, i.e. there is no entanglement in the multipartite system. If fc = 2 the state is called 
biseparable, if it is not biseparable, it is called genuinely multipartite entangled. Genuine multipartite entanglement 
is found to be a key ingredient for many quantum algorithms, see e.g. Ref. [23|. In finite-dimensional systems, there 
exist different inequivalent classes of genuine multipartite entanglement, e.g. the GHZ-class, the M^-class or the 
Dicke-class. Such substructures are only known for very special and rather simple systems (see e.g. [1, [13, HH]) and 
it is not known how this generalises for continuous variable systems. 

Example: A tripartite pure state can be either fully separable, i.e. k ^ n with the 3-partition (xi|x2|a;3), 
or 2-separable (biseparable, i.e. k=2) with one of the three possible partitions (xia;2|a;3), {xiX3\x2) or {xi\x2X3), or 
genuinely multipartite entangled (fc = I) with the partition (a;ia;2a;3). 

For mixed states, we can extend this definition in a straightforward way. A mixed n-partite quantum state 
has the general form 

/oo 
d"a; (Px' p{xi,x[,X2,x'2, ...,Xn,x'^) \xi){x[ \ ® \x2}{x'2 \ «) • • • \xn}{x'J (3) 
-oo 

It is called fc-separable, iff it has a decomposition into fc-separable pure states, i.e. iff it can be written as a convex 
combination of pure fc-separable states: 

p = j dapa \^cc){'^a\ (4) 

where Pa is a probability distribution (i.e. > and da Pa — 1) and \^ a) is fc-separable for all a. 
Note that a fc-separable mixed state may not be separable with respect to any partition, since the pure states 
in its decomposition may be separable with respect to different fc-partitions (i.e. \^a) may split into different 
partitions for different a). The concept of fc-separability is of high impoertance to quantum information theory, 
since many of its applications rely on specific kinds of states (in particular on genuinely multipartite entangled states) . 



III. CRITERION FOR K-SEPARABILITY 

In Refs. [21I, [l^l, a framework for constructing very general separability criteria for finite Hilbert spaces was 
introduced. We now extend this framework such that we can apply it to continuous quantum systems. The main 
result of this paper is an inequality, given in the following theorem, which is satisfied for all fc-separable states, such 
that any violation implies that the state under investigation is not fc-separable. 

Theorem: All fc-separable states p satisfy the inequality 

(k \ik 
l[{'f\plp^'Pc..m] <o (*) 

where |4>) = (g) |02) is an arbitrary fully separable state, Pq^ is a permutation operator permuting the ai-th 
elements of |0i) and \4>2), Ptotai wholly permutes the and the sum runs over all fc-partitions {a}. 
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Proof: Firstly, observe that the left hand side of the inequality (j*J is a convex quantity. Thus, its validity 
for mixed states follows from its validity for pure states. To prove the latter, let us assume w.l.o.g. that a given pure 
state p = \^){^\ is separable with respect to the partition a' . Therefore, the permutation operators corresponding 
to this partition do not change the two copy state 

Furthermore, note that the total permutation acts as 

PtotaMl)®\h) = \h)®\^l) (6) 

then the inequality reads 

Cfe \ A 

Xl{^plp'^^p^^'^)\ <o. (7) 

It follows from the Cauchy-Schwarz inequality that the first two terms cancel (because the pure state p is a projector). 
Since the remaining sum has strictly nonnegative terms with an overall negative sign, the whole inequality has to be 
satisfied for p which are separable with respect to the partition a' . □ 



IV. OPTIMISING THE DETECTION QUALITY 



Evidently, the detection quality of inequality Q strongly depends on the choice of the separable state | $) . Unlike 
in the case of discrete systems, numerical optimisation is quite out of the question, as the number of optimisation 
parameters would be infinite. However, there is a quite intuitive way of choosing |<i>) effectively. Given the state p in 
question, |<i>) should satisfy the following conditions: 

CI: 1$) has to be fully separable, i.e. |$) = 0"^^ (g> {g)"^^ |(/)2i). 

C2: {(pi) and |02) should be orthogonal in each subsystem, i.e. {4>ii\(j>2i) — i- 

C3: Each \(f>ji) should be sharp, i.e. = dx S{x — Xij)\x) — \xij) for some Xij. 

C4: {(pi) and |02) should be chosen such that |(0i|p|02)| is maximal. 

Let us illustrate the background of these conditions. 

Condition CI is necessary for the criterion (ineq. (j*J) to remain valid for all fc-separable states and thus is rather a 
technical requirement. 

Condition C2 guarantees that the permutation operators used in the criterion Q have maximal impact, such 
that a maximal violation of the inequality can be achieved. Each pair \4>ji) {j ~ 1,2) is responsible for detecting 
entanglement in the i-th subsystem. Thus, non-orthonormality of this pair of vectors means non-optimal detection 
of entanglement in this subsystem. 

Condition C3 stems from the fact, that any average is always lower than (or equal to) its highest contribution. Any 
distribution containing more than one element leads to an averaging in the scalar products in Q, which can never 
increase violation of the inequality, but will in general decrease it. 

Condition C4 incorporates the dependance of |$) on the investigated state and thus most of the subtleties involved. 
The first term in ineq. Q (and the only positive one) is the absolute value of the scalar product {(f)i\p\(f)2)- Evidently, 
for the inequality to be violated, this product needs to be as big as possible. This corresponds to chosing |0i) and 
102) as two of the main contributing vectors of the investigated state. 

Although this procedure might not be unique (i.e. might not unambiguously lead to a definite choice of |<&)), it 
reveals the optimal structure of 1$) and thus drastically reduces the number of optimisation parameters to a finite 
set, which can be optimised numerically (e.g. by means of the method introduced in [26| ) or even analytically. This 
will be demonstrated in the next section. 
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V. EXAMPLES 



Consider the family of tripartite states (i.e. n=3) 



iV 

(|xi - A) (g) |x2) (g) |a;3 + A) + |a;i) ® \x2 + A) (g) \x-i - A) + jxi + A) (g) |x2 - A) » jxa) + (8) 
|xi + A) ® 1x2) (g |a;3 - A) + |a;i) ® \x2 - A) ® jsg + A) + \xi - A) <g) \x2 + A) <g) [xg)) 

where 

3 A2(e + 5g) A2(e + 5g) A2(e + 9iT) A^Ce + lStr) ^ g-2A(Ag + 6g) , , ^ 

iV = 27r2eVCT(e ~ + 2e ~ + 2e + 4e ^ + 2e i^^^ + 2e + 2V2e ^ ) (9) 

is a normalisation constant and all parameters a, e and A are larger than or equal to zero. Note that 

lim(27re)-5e~^ = S{x) (10) 

e— fO 

is the Dirac delta function, such that in this limit, for A = the state \uj) is a generalisation of the Gi? Z-state to 
continuous systems: 

1 f°° 

lim |w)a=o = — ^ / dx e \x) (g \x) g) \x) (11) 

-oo 



while for A > it is a generalisation to the W-state, since it is genuinely multipartite entangled and has entangled 
reduced density matrices. 

In the following we are interested in the introduced tripartite states mixed with Gaussian distributed noise 



^2^^2^^2 

Pmix = {2n6y^ I (fxe—^~^^\xi){xi\(S)\x2}{x2\<S)\x3}{x3\ (12) 



i.e. the state under investigation is 

p ^ p \uj) {uj\ + {1 - p) pmix ■ (13) 



A. GHZ-like state (A = 0) 

In the case A = 0, the state \uj) assumes the comparatively simple form 

\uj) = ^= d'^xe'^e ~- |a;i) g) |x2) g) jxa) (14) 

V N J -oo 

As this state's distribution function has its peak at xi = X2 = X3 = and is correlated such that all three coordinates 
are most likely to be very close to each other (or even equal for e — ?> 0), |$) is best chosen to be of the form 

1$) = |0i) g) 102) with 

l^i) = ko) ® |xo) g) |xo) (15) 
102) = I - a;o) g) I - xo) g) I - xo) 

for some xq 7^ (due to condition C2). The optimal value for xq can be obtained from analytical or numerical 
optimisation. Since for cr this state becomes separable (as can be seen directly by substituting eq. (fTO|) in eq. 
p4p ). we are only interested in cr > and thus can, without loss of generality, define cr = 1 and measure all lengths 
in units of ^/a. 

In order to estimate the best value for xo, we first assume a pure state, i.e. p = 1 and investigate the detection 
behaviour of the criterion Q for different values of the two remaining parameters e and xq (as illustrated in Fig. [1} . 
It can be seen that the optimal choice of xq depends on e only slightly and is best chosen in the range 0.7 < xq < 1.2. 
In the further study of this state, we will chose xq = 1, such that only parameters of the state remain. 
Now, we can investigate the mixed state case (i.e. p <1). In Fig. [5] the detection range of our criterion (Q for /c = 2 
and fc = 3 is illustrated . 
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FIG. 1: (Colour online) Illustration of the detection parameter space of the criterion (Q for k — 2 and the state \uj), Eq. (|14|l . 
with p = 1 and A = 0. The state is detected to be genuinely tripartite entangled wherever the graph is above zero (i.e. in the 
green plotted region). It can be seen that the optimal choice for xo is between 0.7 and 1.2, depending only slightly on the value 
of e. 
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FIG. 2: (Colour online) Illustration of the detection quality of ineq. for k — 2 (red) and fc = 3 (blue) for the state (|14|) for 
a — 1, Xo = 1 and different values of p, e and 5. For p = 1 the whole state space is detected to be entangled {k < 2) and to be 
genuinely multipartite entangled (k — 1) for e < 4.648. For lower values of p, still very large portions of the state space are 
detected to be genuinely multipartite entangled (red regions, = 2) or at least entangled (red and blue regions (fc > 2)). 



B. W-like state (A = 1) 

Let us now investigate a more general state, namely the case A 7^ 0. Without loss of generality, we can set A = 1 
and thus measure all lengths in units of A. The optimal choice of |<I>) of course has to be a generalisation of the one 
used in the previous section, which coincides with it for A = 0. The choice is quite straightforward, since the state ^ 
contains six terms. Combination of any two of them leads to the desired high-magnitude off-diagonal density matrix 
element, the only restrictions being conditions Cl-CS. This leads to 

1$) = (8) \(l)2) with 

= \xq + A) (E) \xq) (E) \xo - A) (16) 
|?!'2) = I - Xo - A) (g) I - a:o + A) ® I - Xo) 

In this rather complicated case, numerical optimisation is necessary for achieving optimal detection quality, since 
the optimal xq depends strongly on the other parameters (more than in the previous case with A = 0). However, 
even without numerical optimisation and using only two different choices of xq (namely xq = and xq = 1.5) the 
detection range of our criterion |*]) is very wide, as illustrated in Fig. [3] 



C. Non-Gaussian States 



The above two examples belong to the class of Gaussian states, which have always been the main focus of research 
in the field of continuous variable entanglement (see e.g. Refs. [27H30|). However, recently questions regarding non- 
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FIG. 3: (Colour online) Detection quality of the inequality (Q for the family of the tripartite states ((8| with A = 1, 5 = 1 and 
xo = {0, 1.5} (where a;o = detects better if cr < 1, even for very small p, but ceases to detect for large e; and xo ~ 1.5 detects 
better for a > 1 and still detects a considerable amount of entanglement for large e). The red area is detected to be genuinely 
multipartite entangled (2- inseparable) and the blue and red areas are detected to be entangled (3-inseparable). 



Gaussian entangled states have attracted a lot of interest within the scientific community. Since the criterion presented 
in this paper is - unlike most other criteria for continuous variable entanglement - not tailored specifically for Gaussian 
states, we will also show its detection quality for non-Gaussian states. 
Consider the state 



where 



1 

iV2 



with 



p^p\uj){uj\ + {l-p)pmix (17) 

^ J (fxQ{P - |xi|)e(e - |xi - X2|)e(e - \xi - X3\)\xi)\x2)\x^) (18) 
d^xe(<5 - \xi\)Q{5 - \X2\)Q{5 - |x3|)|xi)(xi| ® |a;2)(a;2| ® \x:i){x^\ (19) 

iVi = 8e^/3 N2 = 8(5^ /3, e, 5 > (20) 

where Q{x) is the Heaviside- function. This state is a non-Gaussian modification of the GHZ-like state ((T4)) . 
Using (jlSp and chosing xq appropriately, the criterion (gj) reads 

{ if § < /3 and ^ < /3 

0< ifi</?<5 (21) 

[ else 

where the first condition is always positive for p > Q and thus indicates genuine multipartite entanglement, the yield 
of the second condition depends on the parameters and the third condition can never be violated. In particular, 
the first condition detects a large portion of the state space to be genuinely multipartite entangled already for 
infinitesimal p > Q. 

Since the state discussed above is only a rather simple example and not very close to experimental realisa- 
tion, let us illustate the detection quality of our criterion for another non-Gaussian state, which is more likely to be 
implemented experimentally [sil [33 ] , namely 

|cj) = aia2a3^= / (Tx e e 2. \xi) ® \x2) ® \xz) (22) 



where the a; are annihilation operators and 

N ^iT^^/^a^'^e" + &ea + l^a^) (23) 



mixed with Gaussian noise, as in (flSj) with (fT2|) . This state also represents a modification of p4|) and can therefore 
be detected by the same choice of 1$). 

Without loss of generality, we set a — 1 and thus measure all lengths in units of y/a. For S < the whole state 
space with p > is detected to be genuinely multipartite entangled, independantly of e. For S > still large areas 
of the state space are detected, as illustrated in Fig. 
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FIG. 4: (Colour online) Illustration of the detection parameter space of the criterion Q for k = 2 and the state p, Eq. (|13p 
where \uj) as in (|22|l . with a = 1 and 5 = 3. The curve indicates the critical proportion p, i.e. the detection threshold for 
genuine multipartite entanglement. 



VI. EXPERIMENTAL IMPLEMENTATION 



Since it is quite important for multipartite entanglement criteria to not only work in theory but also to be imple- 
mentable experimentally (without resorting to a full state tomography, since this would mean an infinite number of 
required measurements), let us illustrate how this can be done for our criterion, ineq. for the pure example state 
from section IV Al the GHZ-like state ([Til) . Generalisation to other states is straightforward, but might be rather 
cumbersome. 

While the criterion detects optimally for sharp 1$), such states cannot be measured physically, since detectors always 
have a finite size. We will thus assume 



with 



\<i>) = \a)\a)\a)\l3m\l3) (24) 



\a)^^Jdxe(^x-xo + ^^e(^xo + ^-x] \x} (25) 
^ dx Q ( X + xo + ^] Q ( + ^ — X \ \x) 



where 8(x) is the Heaviside distribution and ^ is the size of the detector, for example a charge-coupled device (CCD). 
The (2" — 1) scalar products needed for computation of ineq. thus take forms like 



7re 




(26) 



where 



Erf(a) I dx e"^" (27) 



/TT 



is the Gaussian error function. These integrals can easily be computed numerically, given the parameters used in the 
experimental setup, which allows for a simple prediction of measurement outcomes. We now explicitly show how to 
write ineq. ^ in terms of expectation values of local observables in the exemplary three particle case. To that end 
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let us first define the following local observables constructed from the finite sized detectors 

a, = |a)(/3| + |/3)(a| 
ay^i\a){p\-i\P){a\ (28) 
a, = |a)(a|-|/3)(/3|, 

which are the Pauli operators of the two dimensional subspace spanned by \a) and \ jS). Using the short hand notation 

ijk := (o-j (g) cTj (g) crj)p , (29) 
where ai := 1, we can explicitly rewrite ineq.Q for k — 2 and n = 3 as 

— \xxx — yyx — yxy — xyy + i{yyy — xxy — xyx — yxx) \ — 
8 

-{^J{lll + zzl - zlz - Izz + llz - Izl - zll + zzz){lll + zzl - zlz - Izz - llz + Izl + zll - zzz) + (30) 
8 

-/(111 - zzl + zlz - Izz + llz - Izl + zll - zzz)(lll - zzl + zlz ~ Izz - llz + Izl - zll + zzz) + 
\/(lll - zzl - zlz + Izz + llz + Izl - zll - zzz)(lll - zzl - zlz + Izz - llz - Izl + zll + zzz)) < 0. 

It is also possible to decompose the inequalities in terms of local expectation values for larger n or different A: in a 
straightforward way. This however yields rather cumbersome expressions, which is why they are not presented then 
here in full detail. 

Experimental measurement uncertainties can be estimated by means of the Gaussian law of error propagation, which 
states that the measurement uncertainty S of a function / of several arguments Xi is given by 

(31) 

where are the respective measurement uncertainties of the Xi . 

We will assume that all expectation values Xi = {ai\p\ai) underlie the same relative uncertainty (i.e. Q / {ai\p\ai) = C, 
is independant of i), such that only two uncertainty parameters remain, namely Q and o, the latter being the absolute 
uncertainty of the first term in ineq. Q. Now, the measurement uncertainty 5 of the whole inequality is given by 

-'-^'-Y.(if[ '^c) " - i A.f'^e (32, 




where the ineqaHty follows from the fct that a geometric mean is maximal whenever all its factors are equal, and 

/ 1 \k—i-n—l 



is the number of all fc-partitions of an n-partite system. 

In our above examples (i.e. if n = 3 and = 2 or fc = 3), the second term in ([32)) is much smaller than the first, 
such that the measurement uncertainty of the complete expression Q is approximately equal to the uncertainty of 
its first term: S o, which makes this kind of experimental uncertainty easy to deal with. 

Another kind of complication that is to be expected in experimental realisations (e.g. in quantum optics) are 
imperfect detectors. These correspond to nonperfect projective measurements, i.e. each scalar product (ck|p|/?) is 
multiplied by some factor < r < 1. Since the whole inequality is now linear in r, this does not alter the violation 
or nonviolation of the inequality, but reduces the magnitude of violation linearly. 



VII. SUMMARY 



We present a criterion for fc-separability in multipartite continuous variable systems. It is an inequality which 
is satisfied for all fc-separable states, i.e. any violation implies that the state is not ^-separable. The criterion 
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particularly allows to distinguish between biseparable states (which can only be used in few applications) and 
genuinely multipartite entangled ones (which are a basic building piece for several applications of quantum infor- 
mation theory which go beyond the potential of classical systems). We show how the inequality can be optimised 
by chosing an appropriate state |$) (for which we give four explicit conditions) and thus being left with a reduced 
(finite) number of optimisation parameters, which can be computed. We analyse two different families of states, 
which may be considered to be generalisations of the most famous genuinely multipartite states in finite quantum 
information theory, the GH Z-type entangled states and the M^-type entangled states. Our criterion easily detects a 
large parameter space of entangled states when mixed with Gaussian noise. Since no comparable criteria exist, it can 
not be said how tight these detection thresholds are. 

Moreover, we explicitly show how the developed criterion for fc-separability in multipartite continuous variable 
systems can be rewritten by local expectations values, thus how this criterion can be experimentally implemented. 
In summary, we presented a computable criterion for detecting fc-inseparability (and particularly genuine multipartite 
entanglement) in continuous variable systems which can be experimentally implemented by finitely many local 
observables. 
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